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Abstract 



The field theoretic renormalization group and operator product expansion 
are applied to the problem of a passive scalar advected by the Gaussian non- 
solenoidal velocity field with finite correlation time, in the presence of large- 
scale anisotropy. The energy spectrum of the velocity in the inertial range 
has the form E(k) oc k 1 " 6 , and the correlation time at the wavenumber k 
scales as k~ 2+v . It is shown that, depending on the values of the exponents 
e and ij, the model exhibits various types of inertial-range scaling regimes 
with nontrivial anomalous exponents. Explicit asymptotic expressions for the 
structure functions and other correlation functions are obtained; they are rep- 
resented by superpositions of power laws with nonuniversal amplitudes and 
universal (independent of the anisotropy) anomalous exponents, calculated to 
the first order in e and 77 in any space dimension. These anomalous exponents 
are determined by the critical dimensions of tensor composite operators built 
of the scalar gradients, and exhibit a kind of hierarchy related to the degree of 
anisotropy: the less is the rank, the less is the dimension and, consequently, 
the more important is the contribution to the inertial-range behavior. The 
leading terms of the even (odd) structure functions are given by the scalar 
(vector) operators. The anomalous exponents depend explicitly on the degree 
of compressibility. 

PACS numbers: 47.10.+g, 47.27.-i, 05.10.Cc 
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I. INTRODUCTION 



The investigation of intermittency and anomalous scaling in fully developed turbulence 
remains essentially an open theoretical problem. Both the natural and numerical experi- 
ments suggest that the deviation from the predictions of the classical Kolmogorov-Obukhov 
theory |]J2| is even more strongly pronounced for a passively advected scalar field than for 
the velocity field itself; see, e.g., |||J and literature cited therein. At the same time, the 
problem of passive advection appears to be easier tractable theoretically: even simplified 
models describing the advection by a "synthetic" velocity field with a given Gaussian statis- 
tics reproduce many of the anomalous features of genuine turbulent heat or mass transport 
observed in experiments. Therefore, the problem of passive scalar advection, being of prac- 
tical importance in itself, may also be viewed as a starting point in studying anomalous 
scaling in the turbulence on the whole. 

Most progress has been achieved for the so-called rapid-change model, introduced by 
Kraichnan |J: for the first time, the anomalous exponents have been calculated on the 
basis of a microscopic model and within regular perturbation expansions; see, e.g., f|-|§ and 
references therein. 

In the original Kraichnan model, the velocity field is taken to be Gaussian, isotropic, 
incompressible and decorrelated in time. More realistic models should involve anisotropy, 
compressibility and finite correlation time. The recent studies have pointed up some signif- 
icant differences between the zero and finite correlation-time problems [|TU|-|i"2]] and between 
the compressible and incompressible cases [ 13 —pT 7 ? ] . It is noteworthy that the nonsolenoidal 
velocity field remains nontrivial in the one dimensional case, which is more accessible to nu- 
merical simulations and allows interesting comparison between the numerical and analytical 
results; see 



Another important question recently addressed is the effects of large-scale anisotropy 
on inertial-range statistics of passively advected fields [p|,^,|8|, p!2|JT9|] and the velocity itself 



20| . These studies have shown that the anisotropy present at large scales has a strong 



influence on the small-scale statistical properties of the scalar, in disagreement with what 
was expected on the basis of the cascade ideas On the other hand, the exponents 

describing the inertial-range scaling exhibit universality and hierarchy related to the degree 
of anisotropy, which gives some quantitative support to Kolmogorov's hypothesis on the 
restored local isotropy of the inertial-range turbulence |T2] , p!9|j20|] . 

In this paper, we apply the field theoretic renormalization group (RG) and operator 
product expansion (OPE) to the problem of a passive scalar advected by the Gaussian self- 
similar nonsolenoidal velocity field with finite correlation time and in the presence of an 
imposed linear mean gradient. 

Detailed exposition of the RG and OPE approach to statistical models of fully developed 
turbulence and the bibliography can be found in |2~T| , |2"2"|] . The incompressible version of the 
model in question is discussed in [|1^| in detail; below we concentrate on the features specific 
to the nonsolenoidal velocity field. 

The plan of the paper is the following. In Sec. [TT], we describe the model and give its 
field theoretic formulation. In Sec. QI, we discuss the ultraviolet (UV) renormalization of 
the model, derive the RG equations and calculate the basic RG functions (beta functions 
and anomalous dimensions) in the one-loop approximation. The RG equations possess 
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three nontrivial infrared (IR) stable fixed points, which establishes the existence of three 
different IR scaling regimes (Sec. IV). The solution of the RG equations for an equal-time 
two-point correlation function is given, which determines its dependence on the UV scale. 
In Sec. [V], we discuss the renormalization of certain composite operators: powers of the 
scalar field and tensor operators built of the scalar gradients (the former enter into the 
structure functions, and the latter determine the anomalous exponents). In Sec. [VT|, we 
present explicit asymptotic expressions for the structure functions of the scalar field and 
other correlation functions and discuss the hierarchy of the anomalous exponents associated 
with the anisotropic contributions. In Sec. |V1J| , we briefly discuss an alternative model of the 
passive advection by the compressible fluid and present respective inertial-range expressions. 
The results obtained are briefly reviewed in the Conclusion. 



II. DESCRIPTION OF THE MODEL. THE FIELD THEORETIC FORMULATION 

The advection of a passive scalar field in the presence of an imposed linear gradient is 
described by the equation 

d t 6 + (yd)6 = u d 2 6 - (hv). (2.1) 

Here 9{x) = 0(t,x) is the random (fluctuation) part of the total scalar field G(x) = 8(x) + 
(hx), h is a constant vector that determines distinguished direction, h> is the molecular 
diffusivity coefficient, d t = d/dt, di = d/dxi, and d 2 = <9j<9j is the Laplace operator. The 
velocity field v(x) = {vi(x)} obeys a Gaussian distribution with zero mean and correlator 

(Viix^ix')) =J^I 7^p{^i( k ) + «<3i;(k)} D v (u, k) exp[-i(t - f) + ik(x - x')]. 

(2.2) 

Here -Py(k) = 5ij — kikj/k 2 and Qij(k) = kikj/k 2 are the transverse and the longitudinal 
projectors, respectively, k = |k|, and d is the dimensionality of the x space, a > is a free 
parameter. For the function D v we choose 

D v (u;,k) = ^f —j. 2.3 

us 1 + [wo^o k 2 v \ z 

For the energy spectrum we then obtain E{k) ~ k d ~ x j duD v (u, k) ~ g^v^ k 1 ^ 6 . Therefore, 
the coupling constant go > and the exponent e describe the equal-time velocity correlator 
or, equivalently, the energy spectrum, while the constant uq > and the exponent rj are 
related to the frequency uo ~ UqVq k 2 ~ n , characteristic of the mode k. The exponents e and rj 
are the analogs of the RG expansion parameter e = 4 — d in the theory of critical behavior, 
and we shall use the traditional term u e expansion" for the double expansion in the e-rj 
plane around the origin e — r\ — 0, with the additional convention that e = 0(r]). The IR 
regularization is provided by the cut-off in the integral (|2.2| ) from below at k ~ m, where 
i = 1/m is the outer turbulence scale. Dimensionality considerations show that the coupling 
constants go ; u o are related to the characteristic UV momentum scale A by 

g ~ A £ , u ~ A". (2.4) 
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The model contains two special cases that possess some interest on their own. In the 
limit u — > oo, g' Q = go/u = const we arrive at the rapid-change model: 

D v (u;,k)^g' u k- d -^ C = e-ri, (2.5) 

and the limit uq — > 0, go = const corresponds to the case of a "frozen" velocity field: 

D v (u,k) -> g ^k- d+2 ' E 7r5(u), (2.6) 

when the velocity correlator is independent of the time variable t — t' in the t representation. 

The stochastic problem (|2.1| ), (|2.2|) is equivalent to the field theoretic model of the set 
of three fields $ = {9, 9', v} with action functional 

= & [-d t 6 - (yd)9 + u d 2 6 - (hv)] - vD~\/2. (2.7) 

The first four terms in Eq. (|2.7f) represent the Martin-Siggia-Rose-type action (see, e.g., 



Refs. p2|p3| ) for the stochastic problem ( |2.1| ) at fixed v, and the last term represents the 
Gaussian averaging over v. Here D v is the correlator ( |2.2|) , the required integrations over 
x = (t, x) and summations over the vector indices are implied. 

The formulation fl2.7Q means that statistical averages of random quantities in the stochas- 
tic problem fl2.1|) , (|2.2|) coincide with functional averages with the weight expS'($), so that 
generating functional of total and connected [W(A)] Green functions are represented 

by the functional integral 

G(A) = exp W{A) = jv® exp[S($) + A$] (2.8) 

with arbitrary sources A(x) in the linear form 

A$ = J dx[A e {x)6{x) + A e \x)6\x) + AJ(x) Vi (x)}. (2.9) 

The model ( |2.7| ) corresponds to a standard Feynman diagrammatic technique with the 
triple vertex —8'(vd)6 = 9'VjVjO with vertex factor Vj = —ikj, where k is the momentum 
flowing into the vertex via the field 9, and the bare propagators 

{99'), = (9'9)* = (-icu + u k 2 )-\ {99), = {99') h t h 3 {v lVj ) {9'9) , 

{9v i ) = -(^ofcjMo, {9'9') = 0, (2.10) 

where hi is a component of the vector h and the bare propagator {viVj)o is given by Eq. 

(El- 

The magnitude h = |h| can be eliminated from the action (|2.7|) by rescaling of the 
scalar fields: 9 — > h9, 9' — > d'/h. Therefore, any total or connected Green function of the 
form {9(x\) ■ ■ ■ 9(x n ) 9'(y\) ■ ■ •9'{y^)) contains the factor of h n ~ p . The parameter h appears 
in the bare propagators ( [2.10| ) only in the numerators. It then follows that the Green 
functions with n — p < vanish identically. On the contrary, the 1-irreducible function 
(9(xi) ■ ■ ■ 9(x n ) 9'(yi) ■ ■ ■ 9'{y p ))i_ ir contains the factor h p ~ n and therefore vanishes for n—p > 
0; this fact will be relevant in the next Section in the analysis of the renormalizability of the 
model. 
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III. UV RENORMALIZATION. RG FUNCTIONS AND RG EQUATIONS 



The analysis of UV divergences is based on the analysis of canonical dimensions. Dy- 
namical models of the type (|2.7|) , in contrast to static models, have two scales, i.e., the 
canonical ( "engineering" ) dimension of some quantity F (a field or a parameter in the action 
functional) is described by two numbers, the momentum dimension dp and the frequency 
dimension dp. They are determined so that [F] ~ [^"^[T] - ^ where L is the length scale 
and T is the time scale. The dimensions are found from the obvious normalization conditions 
d\ = —d^ = 1, d% = d£ = 0, d k w = d\ = 0, d^ = —df = 1, and from the requirement that 
each term of the action functional be dimensionless (with respect to the momentum and 
frequency dimensions separately). Then, based on d k F and dp, one can introduce the total 
canonical dimension dp = d F + 2dp (in the free theory, dt oc <9 2 ), which plays in the the- 
ory of renormalization of dynamical models the same role as the conventional (momentum) 
dimension in static problems. 

In the action ( |2.7| ), there are fewer terms than fields and parameters, and the canonical 
dimensions are not determined unambiguously. This is of course a manifestation of the fact 
that the "superfluous" parameter h = |h| can be eliminated from the action; see Sec. [TI]. 
After it has been eliminated (or, equivalently, zero canonical dimensions have been assigned 
to it), the definite canonical dimensions can be assigned to the other quantities. They are 
given in Table |, including the dimensions of renormalized parameters, which will appear 
later on. From Table | it follows that the model is logarithmic (the coupling constant go 
is dimensionless) at e — 0, and the UV divergences have the form of the poles in e in the 
Green functions. 

The total canonical dimension of an arbitrary 1-irreducible Green function V = 
($ ■ ■ • $)i_i r is given by the relation 

d v = d r + 2d£ = d+ 2 - N$d$, (3.1) 

where iV$ = {Ng, Ngi, N v } are the numbers of corresponding fields entering into the function 
T, and the summation over all types of the fields is implied. The total dimension dr is 
the formal index of the UV divergence. Superficial UV divergences, whose removal requires 
counterterms, can be present only in those functions T for which dr is a non-negative integer. 
Analysis of the divergences should be based on the following auxiliary considerations: 

(i) All the 1-irreducible Green functions with Ng> < Nq vanish; see Sec. |J. 

(ii) If a number of external momenta occurs as an overall factor in all the diagrams of a 
given Green function, the real index of divergence d' r is smaller than dr by the corresponding 
number of unities (the Green function requires counterterms only if d' T is a non-negative 
integer). In the model fl2.7|) , the field 9 enters into the vertex 9'{yd)9 only in the form of 
the derivative, which decreases the real index of divergence: d' r = dr — Ng. This also means 
that 9 enters into the counterterms only in the form of the derivative 89. 

From the dimensions in Table | we find dr = d + 2 — N v + Ng — id + l)Ngr and d' v = 
d+2 — N v — (d+l)Ngi. Bearing in mind that N$> > Ng we conclude that for any d, superficial 
divergences can exist only in the 1-irreducible functions (#')i_i r with dr = d' r = 1, (0'v)i_i r 
with d r = d' r = 0, (0'0)i-ir with d T = 2, d' T = 1, and (0'0v)!_ ir with d r = 1, d' v = 0. 
The corresponding counterterms necessarily reduce to the forms (hd)9 ! , (hv)9', 9'd 2 9 and 
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9'(vd)9, respectively. The first of these has the form of a total derivative, vanishes after the 
integration over x and therefore gives no contribution to the renormalization action: 



S R (<S>) = 9' \-d t 6 + vZ x d 2 Q - Z 2 (yd)9 - Z 3 (hv)l - vD~ 1 v/2. 



(3.2) 



Here and below the dimensionless parameters g, u, and v are the renormalized analogs of 
the bare parameters, /x is the renormalization mass in the minimal subtraction (MS) scheme, 
which we always use in practical calculations, and Zj = Zi(g,u,a) are the renormalization 
constants. 

The original action ( |2.7|) is invariant with respect to the transformations 9 — > 9 + bx, 
h — ► h — b for any constant vector b. This symmetry is preserved by the renormalization, so 
that the combination Z 2 (vd)9 + Z 3 (hv) must enter into the renormalized action as a whole. 
This implies the exact relation Z 2 = Z 3 . 

The inclusion of the counterterms is reproduced by the multiplicative renormalization of 
the velocity field, v — > Z 2 v, and the parameters go, «o, and in the action functional 



z/ = vZ v , 

The constants in Eqs. (|3.2| ) and ( 

Z v = Zi, 



u = u^ v Z u g = gfi 6 Z g . 
are related as follows: 

Z U = Z X , Zg 



z\zf. 



(3.3) 



(3.4) 



The last two relations in Eq. ( |3.4j ) result from the absence of the renormalization of the 
term with D v in ( |3.2|) . No renormalization of the fields 9, 9' and the parameters m, a, h is 
required, i.e., Z e = 1 and so on. 

In the following, we shall not be interested in the Green functions involving the veloc- 
ity field v, so that we can set A v = in Eq. ( |2.9| ). Then the relation S(9,9',Z 2 v,eo) = 
Sr(9, 9', v, e, /i) (where eo is the complete set of bare parameters, and e is the set of renor- 
malized parameters) for the generating functional in Eq. Q2.8Q yields W(A, eo) = 
Wr(A, e, /i). We use to denote the differential operation /x<9 M for fixed e and operate on 
both sides of this equation with it. This gives the basic RG equation: 



V RG W R (A,e,fx) = 0, 
where T> R q is the operation T>^ expressed in the renormalized variables: 

V RG = Vf, + (3 g d g + (3 u d u - j„V v . 



(3.5) 



(3.6) 

In Eq. ( p.6|) , we have written T> x = xd x for any variable x, and the RG functions (the (3 
functions and the anomalous dimensions 7) are defined as 

1f = V^\yiZ f (3.7) 

for any renormalization constant Zp and 

f3 g =V l ,g = g[-e + 2j 1 -2j 2 ] } f3 u = V„u = u[- V + 7l ]. (3.8) 

The relations between (3 and 7 result from the definitions and the relations (O). 
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Now let us turn to the explicit calculation of the constants Z\^ in the one-loop approx- 
imation. They are determined by the requirement that the 1-irreducible Green functions 
(0'0)i_ir and (6'9\)i_i T be UV finite when expressed in renormalized variables [i.e., have no 
singularities for e, T] — > 0]. The first of these functions satisfies the Dyson equation: 

(0'0)i-ir = -iw + v k 2 - Ee>o(w, k), (3.9) 

where Eg's is the self-energy operator represented by the corresponding 1-irreducible dia- 
grams. In the one-loop approximation it has the form 

5Ve = • (3.10) 



Here and below the solid lines in the diagrams denote the bare propagator (66')o from 
Eq. ( |2.10| ), the end with a slash corresponds to the field 9', and the end without a slash 
corresponds to 9; the dashed lines denote the bare propagator (|2~3| ); the vertices correspond 
to the factor Vi, see Sec. [Q]. The analytic expression for the diagram in ( |3.1(J| ) has the form 

, k) _ fdcu' r dq D v (oj',q)[P ij (q) + aQ ij (q)} 

where q = |q| and D(uj',q) is given by Eq. ( |2.3| ), the factor ki(kj + qj) appears from the 
vertex factors V,. We recall that the integration over q in Eq. ( |3.11|) is restricted from below 
at q ~ m. The integration over uj' yields 

E efe (uj, k) = -ki—± / —— (kj + qj)—— — . , 1V2 , a=r • 3.12 

2 7 (27r) d -io; + z/ (q + k) 2 + wo^o? v 

Equation fl3.12| ) gives the explicit expression for the self-energy operator in the first order 
O(go) of the unrenormalized perturbation theory. Now we need to find Eg/g in the order 
0(g) of the renormalized perturbation theory; therefore we should simply replace g$ — > g/j, £ , 
u — > ufi v , uq — > v in Eq. ( |3.12| ). In the bare term v k 2 = vZ^k 2 in Eq. ( |3.9| ) the 0(g) 
contribution to Z\ should be retained. We know that the divergent part of the diagram is 
independent of u, so that we set uo = in what follows. Since the counterterm is proportional 
to k 2 , we expand the integrand in Eq. ( |3.12|) and neglect all the terms higher than k 2 : 



v / nn 7 f dq q 2 £ [P tj + aQij 

= 0, k) = -fe— J -— (kj + 1j) (Q + k)2 + UQ2(u/ 



" ' dq q 2 - d - £ [ 

(27i) d 1 3 % ' (q + k) 2 + uq 2 (fi/q)v 

gifv f dq q- d - £ [P i:j + aQij] j 2a 



~ hkj 2 J (2n) d l + u(n/q)n Y l+u(fi/q)^ J (3 ' 13) 

where ~ denotes the equality up to an UV finite part, and only the contributions even in ki 
are retained in the integrand. Using the relation 

dqf(q)^ = 6 -ffdqf(q) (3.14) 

we obtain 
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i^Qioyjj = 0, k) — k 



2d 



where we have written 



Jo 



— d — £ 



(2n) d 1 + u (/*/?)" ' 



(rf — 1 + a) J + 2a Ji 



(2vr) d [1 + u (n/q)*] 



(3.15) 
(3.16) 



The renormalized function (#'0v)i_i r in the one-loop approximation is represented as 

(^>i- ir = ^ 2 + /\- (3.17) 
Proceeding as above for (#'#)i_ ir we obtain 



(0'^>i-t 



1 2d 



(3.18) 



with Ji from (|3.16|) . The integrals over q in Eq. ( |3.16j ) are easily performed using the 
expansion in u: 



Jo = E(-«)V 

1=0 



Ji = J2(-u) 1 ^ (i + i) 

«=0 



dq 



-d-e-lr, 
d—e—ltj 



<.in 



-u 



where the parameter m arises from the IR limit in the integral over q, C d = S d /(27i) d and 
S d = 2ir d / 2 /T(d/2) is the surface area of the unit sphere in (i-dimensional space. 

The renormalization constants are found from the requirement that the UV divergences 
cancel out in Eqs. (|3~9|), ( |3.18|) . This determines Z 12 up to UV finite contributions; the 
latter are fixed by the choice of the renormalization scheme. In the MS scheme all the renor- 
malization constants have the form "1 + only poles in e, r\ and their linear combinations," 
which gives the following expressions 

, , gC d a 



lv {-u) l {l + l) 



(3.19) 



Zi 



gC d (d-l + a) gC d a 
1 — <->o H ; — <->i, ^2 



2d 



with C d from Eq. ( p . 1 9[ ) and 



00 f-u) 1 



d 



1=0 



2d 



-u) l {l + l) 
e + It] 



(3.20) 



(3.21) 



The one-loop RG functions can be calculated from the renormalization constants ( |3.20| ) 
using the identity = (3 g d g + (3 u d U) which follows from the definitions ( |3.7| ), (|3.8| ) and the 
fact that Z\ y 2 depend only on the charges g, u. Within our accuracy this identity reduces to 
P„ ~ —tDg ~ V'Du- From Eq. (3.2C) using the relations 



eV g + 7]Vu\ gS = 

1=0 

one obtains: 

gC d 



-u 



1 + u' 



fj 



eV g + r ] V u \gS l = Y J {-u)\l + l)= n 

1=0 \ i ~T~ U ) 



7i 



2d(l+u) 



(d-l + a) 



2a 



72 



gC d a 



2d(l + u 



|2' 



(3.22) 



up to corrections of order g 2 and higher. The beta functions are then obtained from Eqs. 
( g2g ) using the relations (p|). 
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IV. FIXED POINTS AND SCALING REGIMES 



It is well known that possible scaling regimes of a renormalizable model are associated 
with the IR stable fixed points of the corresponding RG equations. The coordinates 
of the fixed points are found from the equations 







(4-1) 



with the beta functions given in Eqs. (3.8). The type of a fixed point is determined by the 
the matrix Q = {Q ik = df3 i /dgk} 1 where $ denotes the full set of the beta functions and g k 
is the full set of charges. For IR stable fixed points the matrix Q is positive, i.e., the real 
parts of all its eigenvalues are positive. 

The analysis of the beta functions ( |3.8| ) reveals five fixed points, which we denote FPI, 
FPU, and so on. In order to find the first two of them, it is convenient to introduce the new 
variables w = 1/u, g' = g/u; the corresponding beta functions have the form 



(3 W = V^w = -Pu/u 2 = w[r] - 71], 
A?' = V^gf = Pg/u- g/3 u /u 2 = g'[rj - e + 71 - 2j 2 ], 



and the anomalous dimensions ( |3.22| ) are written as 



7i 



g'C d 



2d(l + w) 



(d-l + a) 



2aw 



72 



g'w C d a 
2d(l + w) 2 ' 



From Eqs. (|4.2|) and ( [4.3|) we find two fixed points. The first point is trivial, 



FPI : 



= gi = 0; 7* 2 



0. 



The corresponding matrix Q is diagonal with the diagonal elements 

Qi = rj, VL 2 = r] ~ e. 

For the second point we obtain 

2d{e - rj) 



FPU : 



0, giC d 



(d-l + a) 



7i = £ - % 7 2 = 0. 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



The corresponding matrix Q is triangular, d g i(3 w = 0, and its eigenvalues coincide with the 
diagonal elements: 



^i = d w (3 w = rj - 71 = 2rj - e, tt 2 = fy/V = £ 



Now we return to the original variables g, u. From Eqs. 
fixed points. The first one is trivial, 



FPIII : 



u. 



g* 



0; 7i, 



o. 



V- (4.7) 
|3.22| ) we find three more 

(4.8) 



The corresponding matrix Q is diagonal with the elements 
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= -t], Q 2 = s. (4.9) 

For the point FPIV we obtain 

FPIV: ,, = 0, *C<=^; 7l . = 72 .- £ /2=£|ll^. (4.10) 

The corresponding matrix Q is triangular, d g if3 u = 0, and its eigenvalues have the form 

^i = <m = -?7 + 7i, tt 2 = d g {3 g = e. (4.11) 
For the last fixed point FPV both the coordinates ()*■, are finite: 

FPV: 9^, = 2d{e-vl 1+ Mv ~ *) 



l+O (d-l + a)' (d-l + a)(277-e)' 

7i*=*7, 7 2 * = ^-^/2- (4.12) 

A cumbersome but straightforward analysis shows that the positivity condition for Q and 
the condition > (required by the physical meaning of u) are satisfied only in the region 
on the £-77 plane specified by the inequalities e > 0, e(d — 1 — a) < 2r)(d — 1). 

It is noteworthy that the some of the expressions given above are exact, i.e., they have 
no corrections of order e 2 , et], and so on, as a consequence of the exact relations between f3 
and 7 in (|3.7| ) and ( |3.8| ). These are all the results for the trivial points, and 7^ 2 for FPU, 
the relations between Qi and 7^ 2 for FPIV and 7^ 2 for FPV. The expression for 7 2 in Eq. 
( |4.1U| ) is exact only for the incompressible case a = 0. 

It is clear from the Eqs. ( |2.5| ), ( |2.6| ) that the critical regime governed by the point FPU 
corresponds to the rapid-change limit ( |2.5| ) of our model, while the point FPIV corresponds 
to the limit of the frozen velocity field; see Eq. ( |2.6|) . Note that the expression for g' t in 
Eq. Q4.10D coincides with the exact expression obtained in [|l7j] directly for the rapid-change 



model, and is therefore also exact. We then expect that all the critical dimensions at the 
point FPU [FPIV] depend on the only exponent ( = e — 77 [e] that survives in the limit in 
question, and coincide with the corresponding dimensions obtained directly for the models 
( p.5|) [ Q2.6Q 1. This is indeed the case; see Eqs. ( |4.15| ) and (|5.13| ) below. To avoid possible 



misunderstandings we emphasize that the limits go, u — > 00 or u — >• are not supposed in 



the original correlation function (|2.3|) ; the parameters go, Uq are fixed at some finite values. 
The behavior specific to the models ( |2.5| ), ( ^.6[ ) arises asymptotically in the regimes FPU, 
FPIV as a result of the solution of the RG equations, when the "RG flow" approaches the 
corresponding fixed point. More detailed discussion of the physical meaning of the fixed 



points can be found in JT^] for the incompessible case. Triviality of the points FPI and 
FPIII implies the absence of anomalous scaling; we shall not dwell on these regimes in what 
follows. 

In Figure I, we show the regions of stability for the fixed points FPI-FPV in the e-rj 
plane (i.e., the regions for which the eigenvalues of the Q matrix are positive), for some value 
of a > 0. The boundaries of the regions are depicted by thick lines. In the approximation 
linear in e, 77, the regions adjoin each other without overlaps or gaps. When a increases, 
the boundary between FPIV and FPV (i.e., the ray 2i](d — 1) = e(d — 1 — a), e > 0) moves 
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upwards and for a — > oo coincides with the boundary between FPIII and FPIV (e = 0, 
f] > 0), so that the region of stability for FPIV shrinks to zero size. Moreover, the fixed 
point ( fQU|) for a — > oo disappears. In order to treat this limit accurately, it is necessary to 



change to the new variable g" = ga and keep it finite as a — > oo. Then the beta function 
(3 g ii = a/3g'> for a — > oo, u = becomes linear, (3 g n = —eg", and has no fixed points other 
than g" = 0. This fact holds valid for all orders of the perturbation theory, which becomes 
clear from the comparison with the models of the random walks in random environment 
with long-range correlations; see Refs. [f24|-|26|| and references therein.^] 



When a decreases, this boundary moves clockwise and for a = coincides with the ray 
e = 2r] > 0, the boundary between FPU and FPIV. Hence, the region of stability for FPV 
turns into the ray e = 2r\ > 0, in agreement with the result obtained in directly for the 
incompressible case. In the latter, the coordinates of the fixed points and the anomalous 
exponents are nonuniversal in the sense that they can depend not only on the exponent 
£ = 2r], but also on the amplitudes g , u ] see [|TU| , |T!2|| . This nonuniversality can be viewed 
as a consequence of the ambiguity of the values of g*, in Eq. (f4.12|) in the limit a — > 0, 
e — » 2r] (the limit of the combination g*Cd/(l + «*) = 2dr}/(d — 1) is unambiguous and 
coincides up to the notation with the result obtained in [12[ directly for the incompressible 
case) . 

Let F(r) be some equal-time two-point quantity, for example, the pair correlation func- 
tion of the primary fields 9, 8' or some composite operators. We assume that F(r) is 
multiplicatively renormalizable, i.e., F = ZpF R with certain renormalization constant Zp. 
The existence of nontrivial IR stable fixed point implies that in the IR asymptotic region 
Ar >> 1 and any fixed mr the function F(r) takes on the form 

F(r) ~ vf A dF (Ar)- AF £(mr), (4.13) 

where d F and dp are the frequency and total canonical dimensions of F, respectively, and 
£ is some function whose explicit form is not determined by the RG equation itself. The 
critical dimension is given by the expression 

A[F] = A F = d F + A^dp + 7^1, (4.14) 

where 7^ is the value of the anomalous dimension ( p.7| ) at the fixed point and A w = 2 — 7* 
is the critical dimension of frequency. From Eqs. (|3.4|) and (|3.7|) it follows 7^ = 71, so that 
for the nontrivial fixed points we obtain 



C for FPU, 

+ 0(e 2 ) for FPIV, (4.15) 



e(d — 1 — a) i 2 



2(d-l) 

7] for FPV 

(we recall that ( = e — 77, see ( p.5| )). The results for FPU and FPV are exact and coincide 
with their analogs for the incompressible case. 



1 Strictly speaking, the model (|2.1| ) in the frozen limit ( |2.6D differs from the model of random 



random walks studied in [24-26], but their basic RG functions coincide; see Sec. VII 



11 



Note that the dimensions d F , dp and A^ of the pair correlator F(r) = (Fi(x)F2(x')) 
are equal to the sums of the corresponding dimensions of the quantities The critical 
dimensions of the fields 9, 9' in our model are found exactly: A g = —1 A e i = d + 1, and for 
the IR scale we have A m = 1 (we recall that all these quantities in the model (|2.7f) are not 
renormalized and therefore their anomalous dimensions vanish identically, jf = 0). Note 
also that the dimensions Ag^ are independent of a and coincide with their analogs for the 
incompressible case. 



V. CRITICAL DIMENSIONS OF COMPOSITE OPERATORS 

Any monomial or polynomial constructed of primary fields and their derivatives at a 
single spacetime point x = (t, x) is termed a composite operator. Coincidence of the field 
arguments in Green functions containing a composite operator F gives rise to additional 
UV divergences, which are removed by the additional renormalization procedure. As a rule, 
the renormalization of composite operators involves mixing, i.e., an UV finite renormalized 
operator is a linear combination of unrenormalized operators, and vice versa. The detailed 
discussion of the renormalization of composite operators in turbulence models can be found 
in pi| , |2"2"|| ; below we confine ourselves to only the necessary information. 

Let F = {F a } be a closed set, all of whose monomials mix only with each other in 
renormalization. The renormalization matrix Zp = {Zap} and the matrix of anomalous 
dimensions = {70/3} for this set are given by 

F a = Y, z ^Ff, lF = Zp'V^Zp, (5.1) 



and the corresponding matrix of critical dimensions Ap = {A Q/ g} is given by Eq. ( (4.14| ), 
in which dp and d F are understood as the diagonal matrices of canonical dimensions of the 
operators (with the diagonal elements equal to sums of corresponding dimensions of all fields 
and derivatives entering into F) and 7^ = 7f(5'*, m *) is the matrix ( |5TT|) at the fixed point 
in question. Critical dimensions of the set F = {F a } are given by the eigenvalues of the 
matrix A^. Owing to the renormalization, the critical dimension associated with certain 
operator F is not in general equal to the simple sum of critical dimensions of the fields and 
derivatives constituting F. 

The simplest composite operators 9 n (x) in the model are UV finite, 9 n = Z n [9 n ] R with 
Z n = 1. It then follows that the critical dimension of 9 n is simply given by the expression 
( |4.14j ) with no correction from 7^ and is therefore reduced to the sum of the critical dimen- 
sions of the factors: A[9 n ] = nA[9] = —n. The proof is almost identical to the analogous 
proof for the Kraichnan model, given in ||, and will not be discussed here. 

An important role will also be played by the operators of the form 

F[n,p} = d n 9---d tp 9(d l 9d l 9) l 1 (5.2) 

where p is the number of the free vector indices and n = p + 21 is the total number of the 
fields 9 entering into the operator; the vector indices of the symbol F[n,p] are omitted. The 
analysis similar to that given in Ref. ]12| shows that these operators mix only with each 



other in renormalization, the corresponding matrices in Eq. ( |5.1| ) are triangular, and the 
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critical dimensions associated with the family ( |5.2| ) are given by A[n,p] = 7*[n,p], where 
7[n,p] is the diagonal element of jp. The "basis" operator that possesses definite critical 
dimension A[n,p] is a p-th rank tensor (traceless for p > 2), given by a linear combination 
of the monomials F[n',p'] with n' < n, p' < p; the coefficients involve the vector h and 
Kronecker delta symbols. 

Now let us turn to the one-loop calculation of the diagonal element Z[n,p] of the matrix 
Zp in the MS scheme. Let T(x; 9) be the generating functional of the 1-irreducible Green 
functions with one composite operator F[n,p] and any number of fields 9. Here x = (t, x) is 
the argument of the operator and 8(x) is the functional argument, the "classical analog" of 
the random field 9. We are interested in the n-th term of the expansion of T(x; 9) in 9(x), 
which we denote T n (x; 9); it has the form 

T n (x; 9) = J dxi - ■ ■ J dx n 9(xi) ■ ■ -9(x n ) (F[n,p](x)9(x 1 ) ■ ■ ■ 6>(x n ))i_ ir . (5.3) 
In the one-loop approximation the functional (p.3[) is represented diagramatically as follows: 




T n = F[n,p} + ^ (5.4) 

The first term is the "tree" approximation, and the black circle with two attached lines in 
the diagram denotes the variational derivative V(x; Xi,x 2 ) = 5 2 F[n,p]/59(xi)59(x2). It is 
convenient to represent it in the form 

8 2 

V(x; xi,x 2 ) = di5(x - xi) dj8(x - x 2 ) n [a h ■■■a i (a 2 ) 1 }, (5.5) 

OdiOaj 

where a« is a constant vector, which after the differentiation is substituted with di9{x). 

The vertex ( |5.5| ) contains (n — 2) factors of d9. Two remaining "tails" 9 are attached to 
the vertices 9'(\d)9 of the diagram ( |5.4j ). It follows from the explicit form of the vertices 
that these two fields 9 are isolated from the diagram in the form of the external factor 
8989. In other words, two external momenta, corresponding to these fields 9, occur as an 
overall factor in the diagram, and the UV divergence of the latter is logarithmic rather than 
quadratic. Therefore, we can set all the external momenta in the integrand equal to zero, 
and the UV divergent part of the diagram (|5.4j) takes on the form 



8 2 

a p a s ^ da • • • a ip (a 2 ) 1 } T ij;PS , (5.6) 



where we have denoted 



lhps J 2n J (2n) d ^ V lH> u 2 + v 2 q 4 ' V ' 

with D v from Eq. (|2~3|). In Eq. (p.7|), we have to substitute g — > gfi £ , u ufi v and 
uq — > v] see the comments below Eq. (|3.11|) in Sec. |T|. We perform the integration over u 
and use the relations ( 3.14|) and 



, £/ \Q.iQ.jQ.i% ^ijSip -\- SipSij + 5n5pj . 
dqf{q) ^ 4 = d{d + 2) ' dq ^" 
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This gives 



T 



ij.ps 



gp e Jo 
2d(d + 2) 



(d+l + a)SijS ps + (a - l)(5 is 5 jp + 8 ip 5 3S , 



(5.8) 



with the integral Jq defined in Eq. ( |3.16 ). 

Substituting Eq. ( ]5.8| ) into Eq. ( |5.6| ) gives the desired expression for the divergent part of 
the diagram (|5.4p . In this expression we have to take into account all the terms proportional 
to the operator F[n,p] and neglect all the other terms, namely, the terms containing the 
factors of 5^ etc. The latter determine non-diagonal elements of the matrix Zp, which we 
are not interested in here. Finally we obtain 



r ~ 



F[n,p] 



gp £ J Q[n,p] 
Ad(d + 2) 



+ 



(5.9) 



where we have written 



Q[n,p]=2p(p-l)(l 



a) 



[n 



p) (d - 1 + a)(n +p + d) + 2a(n + p - 2) (5.10) 



and Jo is defined in Eq. ( 3.16|) . The dots in Eq. ( |5.9| ) stand for the 0(g 2 ) terms and the 
structures different from F[n,p], ~ denotes the equality up to UV finite parts; we also recall 
that n = p + 21. 

The constant Z[n,p] is found from the requirement that the renormalized analog = 
Z~ l \n,p]T n of the function ( |5.9| ) be UV finite (mind the minus sign in the exponent); along 
with the representation (|3.19|) for the integral J and the MS scheme this gives 



Z[n,p] 



Q[n,p] 2 

9Cd AdJdTY) So + 0{9) 



(5.11) 



with Cd from Eq. ( p,19|) and So from Eq. (|3.21|) . For the the anomalous dimension 7 [n, p] 
T)^ In Z[n,p] it then follows (see the comments below Eq. ( |3.21| ) in Sec. fTTl] ) 



l[n,p] 



gC d Q[n,p] 



(u + 1) Ad(d + 2) 
For the nontrivial fixed points discussed in Sec. |TT^ we then obtain 

1 + cO 



(5.12) 



A[n,p] 



Q\np] \ 2C ^ d - 



for FPU, 
for FPIV, 



4(d + 2) [ 2 (/(d-l + a) forFPV 



(5.13) 



up to corrections of order e 2 and so on. The expression ( 5.13|) illustrates the general fact that 
the critical dimensions in the rapid-change and frozen regimes depend only on the exponents 
( and e, respectively. The coincidence of the results for the points FPIV and FPV seems 
to be an artifact of the one-loop approximation. For a = 0, the results of Ref. [|12[ 
recovered. We also note that the dimensions Q5.13 ) have well-defined limit d — > 1 for p 
and 1, i.e., for the scalar and vector operators, and this limit does not depend on a. 



are 
= 
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The result A [2,0] = — ( for the rapid-change regime is in fact exact. This can be 
demonstrated using the Schwinger equation of the formf] 



J V<$>5 [6(x) exp S R ($) + A$] /56'(x) = (5.14) 



In Eq. ( |5.14|) , Sr is the renormalized action ( |3.2|) , and the notation introduced in Eqs. 
Q2.9| ) is used. Equation ( J5.14 ) can be rewritten in the form 



-vZ x F% 0] + ^Z 2 (d iVl )6 2 - Z 2 (hv)9 + ^{-d t 9 2 - Z 2 d t [v t 9 2 ] + vZ x &P} S 

= -A e ,5W R (A)/5A e . (5.15) 



Here (. . .)a denotes the averaging with the weight exp [£#($) + A&], Wr is determined by 
Eq. ( |2.8|) with the replacement S — > Sr, the argument x common to all the quantities in 
( 15.151 ) is omitted, and the notation ( |5.2|) is used. 

The quantity (F)a is the generating functional of the correlation functions with one 
operator F and any number of the primary fields <3>, therefore the UV finiteness of the 
functional (F)a is equivalent to the finiteness of the operator F itself. The quantity on the 
right hand side of Eq. (|5.15|) is UV finite (a derivative of the renormalized functional with 
respect to a finite argument), and so is the operator on the left hand side. We are going to 
prove that the operator F[2, 0] is UV finite on its own. None of the operators other than 
F[2, 0] can admix to F[2, 0] in renormalization (the family (|5.2|) is closed with respect to 
the renormalization, see above). Our operator F[2,0] does not admix to (hv)8 (no needed 
diagrams can be constructed), and to the operators in the curly brackets (they have the 
form of total derivatives, and F[2, 0] does not reduce to this form). The admixture of F[2, 0] 
to the remaining operator (diVi)8 2 is determined by the 1-irreducible function 

((d l v i )9 2 (x)9(x l )9(x 2 )) l ^ = • • +•■■ (5.16) 
which in the one-loop approximation is written as 

- * / W) I Wr^ + ^ U+ffiXd. ■ (5 ' 17) 

see Eqs. (|2.2|) , (|2.10| ). In the rapid-change limit, the function D v (u',q) = D v (q) does not 
depend on the frequency, see Eq. (|2.5j) , and the integration over to' gives the result 

ah f dq 

-qiD v (q), 



(2tt)« 

which vanishes owing to the symmetry (it is important here that the integrand is indepen- 
dent on k). Furthermore, the multiloop diagrams of the function ( |5.16| ) contain effectively 
closed circuits of retarded propagators (99') and therefore vanish (it is crucial here that the 



2 In the general sense of the term, Schwinger equations are any relations stating that any functional 
integral of a total variational derivative vanishes; see, e.g., [22, 23 1 . 
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propagator ( |2.2| ) in the rapid-change case is proportional to the S function in time). There- 
fore, the operator F[2, 0] does not admix to (<9jt>j)# 2 ; it is completely independent, and it 
must be UV finite separately. 

Since the operator uZiF[2,0] is UV finite, it coincides with its finite part, i.e., 
uZiF[2,0] = uF R [2,0], which along with the relation F[2, 0] = Z[2, 0] F R [2, 0] gives 
Z[2,0] = Z-f 1 and therefore j[2, 0] = —71. For the critical dimension we then obtain 
A[2, 0] = 7* [2, 0] = -C exactly (we recall that 7? = 7* = e - rj = £ see Eq. ([Op). Note 
that the UV finiteness of the remaining combination (diVi)8 2 /2 — (hv)6 implies Z 2 — 1 for 
the rapid-change regime (g/u = const, u — > 00), which can be directly verified using the 
explicit expression ( |3.20| ). This is a consequence of the fact that all the nontrivial diagrams 
of the 1-irreducible function (##'v)i_i r also contain closed circuits of retarded propagators 
and vanish, cf. Jl7| . 



For the regimes FPIV and FPV, the integral (|5.17|) and the higher-order diagrams are 
nontrivial, F[2, 0] admixes to (diVi)6 2 , and the proof does not work. Exact results for A [2, 0] 
are then available only if a = 0, when the operator (diVi)8 2 vanishes identically; see |12| . 



VI. OPERATOR PRODUCT EXPANSION AND THE ANOMALOUS SCALING 

The representation ( [4.13| ) for any scaling function £(mr) describes the behavior of the 



Green function for Ar >> 1 and any fixed value of mr. The inertial range corresponds to the 
additional condition that mr << 1. The form of the function £(mr) is not determined by 
the RG equations themselves; in the theory of critical phenomena, its behavior for mr — > 
is studied using the well-known Wilson operator product expansion; see, e.g., This 
technique is also applicable to the theory of turbulence; see plp2| . 



According to the OPE, the equal-time product Fi{x)F2(x') of two renormalized operators 
at x = (x + x')/2 = const and r = x — x' — ► has the representation 

F 1 (x)F 2 (x') =£C a (r)F a (x,t), (6.1) 

a 

where the functions C a are the Wilson coefficients regular in m 2 and F a are all possible 
renormalized local composite operators allowed by symmetry, with definite critical dimen- 
sions A a . The renormalized correlator (Fi(x)F 2 {x')) is obtained by averaging Eq. ( |6 . 1| ) 
with the weight exp Sr, the quantities (F a ) appear on the right hand side. Their asymp- 
totic behavior for m — >• is found from the corresponding RG equations and has the form 
(F a ) oc m Aa . From the operator product expansion ( |6.1| ) we therefore find the following 
expression for the scaling function ((mr) in the representation ( |4.13| ) for the correlator 
(F 1 (x)F 2 (x')): 

(( mr ) = A a (mr) Aa , (6.2) 

a 

where the coefficients A a = A a (mr) are regular in (mr) 2 . 

The quantities of interest are, in particular, the equal-time structure functions 

S n (r) = ([9(t,x)-e{t,x!)] n ), r = |x-x'|. (6.3) 
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For these, the representation ( |4.13| ) is valid with the dimensions d% = and 4f = Ap = 
nAg = —n. In general, the operators entering into the OPE are those which appear in 
the corresponding Taylor expansions, and also all possible operators that admix to them in 
renormalization. The leading term of the Taylor expansion for the function (|6.3|) is given 
by the n-th rank tensor F[n, n] from Eq. ( |5.2[ ). The decomposition of F[n, n] in irreducible 
tensors gives rise to the operators F[n,p] with all possible values of p < n; the admixture of 
junior operators gives rise to all the monomials F[k,p] with k < n and all possible p. Hence, 
the desired asymptotic expression for the structure function has the form 

S n (r) ~ (hr) n E [ C *p H A[M +•••], (6-4) 

k=0P=Pk 



with the dimensions A[k, p] from Eq. (|5.13 ). Here and below pk denotes the minimal possible 
value of p for given k, i.e., pk = for k even and p^ = 1 for k odd; C^ v are some numerical 
coefficients dependent on e, rj, d, a, and on the angle •& between the vectors h and r. The 
dots in Eq. (|6.4j) stand for the contributions of order (mr) 2+0 ' e ' and higher, which arise 
from the senior operators like d 2 9d 2 9. The operators F[k,p] with k > n (whose contributions 
would be more important) do not appear in Eq. fl6.4j) , because they are absent in the Taylor 
expansion of S n and do not admix in renormalization to the terms of the Taylor expansion. 

The mean value of the operator with the dimension A[&,p] is a traceless (for p > 2) 
tensor built of the vector h and Kronecker delta symbols; the contraction with the vector 
indices of the corresponding coefficient C a (r) in Eq. (|6.1|) gives rise to p-th order Legendre 
polynomial P p (cos$), so that the expansion Q6.4Q is analogous to the Legendre polynomial 



decomposition of Ref. [TjJ, and to the decomposition in irreducible representations of the 
rotation group, used in Ref. [pOfl . 

The straightforward analysis of the explicit one-loop expression ( |5.13| ) shows that for fixed 
n, any a > 0, and any nontrivial fixed point, the dimension A[n,p] decreases monotonically 
with p and reaches its minimum for the minimal possible value of p = p n , i.e., p = if n is 
even and p — 1 if n is odd: 

A[n,p] > A[n,p'] if p>p'. (6.5) 

A similar inequalities are satisfied by the critical dimensions of certain tensor operators in 
the stirred Navier-Stokes turbulence, see Ref. [p7| and Sec. 2.3 of f2~2|| , and in the model 
describing passive advection of the magnetic field by the rapid-change velocity in the presence 
of a constant background field ]nj. Furthermore, this minimal value A[n,p n ] is negative 



and decreases monotonically as n increases: 

> A[2k, 0] > A[2k + 1, 1] > A[2k + 2, 0]. (6.6) 

The inequalities ( |6.5| ), ( |6.6|) show that the contrbutions of the tensor operators ( |5.2| ) into 
the asymptotic expression ( |6.4j) exhibit a kind of hierarchy: the less is the rank, the more 
important is the contribution; cf. |12| for the purely incompressible case. 

The leading term of the expression ( |6.4j) for the even (odd) function S n is determined by 
the scalar (vector) composite operator consisting of n factors 89 and has the form 

Sjr) oc {hr) n {mr) A ^. (6.7) 
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It is easily seen from the explicit expression ( 5.13|) that Eq. (|6.7|) for the rapid-change regime 



and even n is in agreement with the results obtained in Refs. fl5| , [16ll for a "tracer" and earlier 
in |§ for d = 1.0 



Expressions similar to Eqs. (|6.4j ), ( |6.7|) can be written down for other correlation func- 
tions, provided their canonical and critical dimensions are known; in particular, the analog 
of the expression ( |6.7|) for the equal-time pair correlator of the operators ( |5.2|) has the form 



(F[n,p\ F[n',p'}) ~ h n+n ' (Ar)- A ^~ A ^(mr) A t n+n '' p "+"' ] , (6.8) 
with the dimensions A[n,p] from Eq. ( |5.13| ); cf. [O for the incompressible case. 

VII. PASSIVE ADVECTION OF A DENSITY FIELD 

So far, we have studied the model ( |2.1|) with the nonlinearity of the form (yd)9, which 
describes a "tracer" in the terminology of Refs. |~5],[16|. Another possibility is to choose 
d(yd), so that the dynamical equation 

d t 9 + 0(v0) = u d 2 6 - 9i(vi(hv)) (7.1) 

for h = becomes a conservation law for 6; see ]r5|JT7[ ]. The "frozen" limit fl2.(j| ) of Eq. ( [7.1| ) 
with h = leads to the model of the random walks in random environment with long-range 



correlations; see [ 24, 26 



For h = 0, the action functionals corresponding to the models ( |2.1| ) and ( |7.1| ) are con- 
nected by the transformation 0^6' and t — ► —t. Since the field h does not enter into 
the diagrams needed for the calculation of the constants Z 12 (see Sec. [HID , t ne basic RG 
functions, dimensions A u g gi, coordinates of the fixed points and their regions of stability 
for the models ( [2.1| ) and ( [T.l|) coincide identically. Note that the expression ( 4.15|) for FPIV 



is in agreement with the result obtained in the model of random-random walks f24|-p6| (up 



to notation and a misprint in Eq. (4.48) of Ref. |25]: 2u in the first line should be replaced 
with (2z/) _1 ; in our notation, A w = 1/u). 

The symmetry 9 — > 6 + const , specific to the model (|2.1|) , ceases to hold for the model 
(7TT), and the real index of divergence for this case has the form d' T = d-p — Ngi (the derivative 



d at the vertex 6'd{y6) can be moved onto the field 6' using the integration by parts). As a 
result, the analysis of composite operators differs essentially in the two cases: in particular, 
the operators 9 n in the model ( |7.1|) require nontrivial renormalization, 9 n = Z n [6 n ] R with 
Z n ^ 1; cf. WJ§ for the rapid-change case. The one- loop calculation of Z n is similar to the 
calculation of the constants Z 12 in Sec. |T| and Z[n,p] in Sec. [V], so that we present only 



the result: 

Z n = l + y - 4 l 1 S + 0{g 2 ), (7.2) 

with Cd from Eq. (|3.19| ) and Sq from Eq. ( p.21|) . Hence, the operators 9 n acquire nontrivial 
anomalous dimensions, 7„ = V^liaZ n = —an{n — l)gCd/^{u + 1) (they should not be 



3 In the notation of Refs. 15,1611, p = a/(d — 1 + a). 
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confused with 7^2 from Sees. |I| |TV|) , and for the corresponding critical dimensions A n = 
—n + 7* we obtain 

cm(w- l)rf f 2C/(d - 1 + a) for FPU and FPV, . , 

An _ H 4 X \e/(d-l) for FPIV, {1 - 6) 

cf. Eq. ( |5.13| ). The ( 2 contribution to 7* and the exact expression for 72 in the rapid-change 
regime are given in Ref . |17] . 



Consider now the equal-time pair correlator {9 p (t, x.)9 k (t, x')). Substituting the rela- 
tions d% = and dp = —(p + k) into the general expression ( 4.13|) gives (9 P 9 ) 



\ ( p+fe )(Ar) Ap+Afc ^ P) fc(mr), with A n from Eq. ( |7.3|) (here and below, we do not display 
the obvious dependence on h). The small mr behavior of the scaling functions is found 



from Eq. (|6.2| ). In contrast to the examples from Sec. |VT|, the composite operators in the 
expansion ( |6.1| ) can involve the field 9 without derivatives. The leading term in Eq. ( |6.2| ) is 
then determined by the monomial 9 p+k , which gives 

(9 p 9 k ) ~ A-( p+k \Ar)- A p- Ak (mr) A i>+K (7.4) 

Note that Eq. ( [7.4| ) contains nontrivial dependence on both the IR and UV scales. 

Now let us turn to the structure functions (|6.3|) in the inertial range Ar >> 1, mr << 1. 
From the expression (|7.4J) it follows 



(7.5) 



S n (r)~A~ n (m/A) A " 1+ £ C kp (Ai 



k+p=n 
k,p^0 



where the coefficients C^p are independent of the scales A, m and the separation r. It is 
obvious from the inequality + A p > A k + P , satisfied by the dimensions ( [7.3|) , that all the 
contributions in the sum in Eq. (|7.5|) vanish in the limit Ar — > 00, so that the leading term 
of the structure function does not depend on r and is simply given by S n ~ (6* n ). 



The situation changes for the purely solenoidal velocity field, a = in Eq. (2.2). In 



this case, the models (2.1) and ( |7.1|) become identical, the dimensions (|7.3| ) become linear 



in n, A n = —n, and the expression ( |7.4| ) becomes independent of r. The operators 9 n 
cannot appear in the OPE for the structure functions. This means that the contributions 
of the operators 9 n to the pair correlators ( fT^ ) cancel out in the functions (|6.3|), and the 
IR behavior of the latter is dominated by the operators constructed solely of the scalar 
gradients. The cancellation becomes possible owing to the fact that all the terms in the 
curly brackets in Eq. ( |7.5| ) become independent of Ar. In this case, the anomalous scaling 
of S n is determined by the critical dimensions of the operators (|5.2| ) from Sec. 0. 



VIII. CONCLUSION 

We have applied the field theoretic renormalization group and operator product expan- 
sion to the model of a passive scalar advected by the Gaussian self-similar nonsolenoidal 
velocity field with finite correlation time and in the presence of large-scale anisotropy, in- 
duced by the linear mean gradient. The energy spectrum of the velocity in the inertial range 
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has the form E(k) oc gok 1-6 , and the frequency at the wavenumber k scales as uok 2 ~ v . It 
is shown that, depending on the values of the exponents e and rj, the model exhibits three 
types of the inertial-range scaling regimes with nontrivial anomalous exponents. 

The explicit asymptotic expressions for the structure functions and other correlation 
functions are obtained; the anomalous exponents, determined by the critical dimensions 
of certain composite operators, are calculated to the first order in e and rj in any space 
dimension. For the first scaling regime the exponents are the same as in the rapid-change 
limit of the model; for the second they are the same as in the model with time-independent 
(frozen) velocity field. For all these regimes, the anomalous exponents are nonuniversal 
through their dependence on a, the degree of compressibility. However, they are independent 
of the amplitudes go, uq, in contradistinction with the case of a local turnover exponent 
(e = 2rj) in an incompressible flow fllO| , |T2"| . 

The most interesting result of the RG analysis is that the anomalous exponents, related 
to the anisotropic contributions to the inertial-range behavior (or, in other words, critical 
dimensions of the tensor composite operators in the corresponding OPE) exhibit a kind of 
hierarchy related to the degree of anisotropy: the less is the rank, the less is the dimension 
and, consequently, the more important is the contribution; cf. for the purely incompress- 
ible case. The leading term for the even (odd) correlation functions is then determined by 
the scalar (vector) composite operators (i.e., those having minimal possible rank). A sim- 
ilar hierarchy is demonstrated by the critical dimensions of certain tensor operators in the 
stirred Navier-Stokes turbulence, see Ref. [p7| and Sec. 2.3 of p2f , and in the model of the 
magnetic field advected passively by the rapid-change velocity It can be viewed as an 
analytic argument in favor of the Kolmogorov hypothesis on the restored local (small-scale) 
isotropy; see the discussion in Ref. [[HJ . 



To avoid possible misunderstandings, it should be emphasized that the large-scale 
anisotropy persists, through the dependence on the imposed field h, for all ranges of mo- 
menta (including inertial and dissipative ranges). Indeed, the correlation function involving 
n fields 9 is proportional to h n , which is an obvious consequence of the linearity of the 
equation (2.1) in 9 and h (see also the discussion in the end of Sec. |l|). It then follows 
that the dimensionless ratios of the structure functions are strictly independent on h and 
can have finite limit as h — > 0, in spite of the fact that the structure functions themselves 
vanish for h = 0. It is noteworthy that all these statements equally hold for any velocity 
field (not necessarily Gaussian or synthetic), provided its statistics is independent of h, in 
agreement with the recent findings of Refs. |||| , where the (derivative) skewness factor has 
been analyzed. 

The explicit inertial-range expressions ( p.4[ ), ( |6.8| ), ( |7.5[ ) show, however, that only the 
amplitudes in power laws are affected by the anisotropy (or by the statistics of the stirring 
force in the original Kraichnan model, or by the initial and boundary conditions in a finite- 
size problem). The exponents are independent of h and can be calculated directly in the 
homogeneous model with h = 0. It follows from the derivation of Eqs. ( |6.4|) , ( |678D , ( [7.51) , 
that this statement holds valid not only for the leading exponents, related to the dimensions 
of scalar operators (|5.2|) with p = 0, but also for all the anisotropic corrections, related to 
the tensor operators with p ^ 0! 

The picture outlined above and in Refs. fi2| , |T9H for passively advected fields (a superpo- 
sition of power laws with universal exponents and nonuniversal amplitudes) seems rather 
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general, being compatible with that established recently in the field of Navier-Stokes tur- 
bulence, on the basis of numerical simulations of channel flows and experiments in the 
atmospheric surface layer, see Refs. |B0 and references therein. In those papers, the veloc- 
ity structure functions were decomposed into the irreducible representations of the rotation 
group. It was shown that in each sector of the decomposition, scaling behavior can be 
found with apparently universal exponents. The amplitudes of the various contributions 
are nonuniversal, through the dependence on the position in the flow, the local degree of 



anisotropy and inhomogeneity, and so on [20 



However, the anomalous exponents of the passive scalar become nonuniversal and acquire 
the dependence on the anisotropy parameters if the velocity field remains anisotropic at small 
scales [p8| . 



ACKNOWLEDGMENTS 

I have benefited from discussions with L. Ts. Adzhemyan, Juha Honkonen and 
M. Yu. Nalimov. I also thank Andrea Mazzino for useful comments on the paper subject. 

The work was supported by the Grant Center for Natural Sciences of the Russian State 
Committee for Higher Education (Grant No. 97-0-14.1-30) and by the Russian Foundation 
for Fundamental Research (Grant No. 99-02-16783). 



21 




22 



REFERENCES 



[1] A. S. Monin and A. M. Yaglom, Statistical Fluid Mechanics (MIT Press, Cambridge, 
MA, 1975), Vol.2. 

[2] U. Frisch, Turbulence: The Legacy of A. N. Kologorov (Cambridge University Press, 
1995). 

[3] K. R. Sreenivasan, Proc. Roy. Soc. London A 434, 165 (1991); K. R. Sreenivasan and 

R. A. Antonia, Annu. Rev. Fluid. Mech. 29, 435 (1997). 
[4] M. Holzer and E. D. Siggia, Phys. Fluids 6, 1820 (1994); A. Pumir, Phys. Fluids 6, 

2118 (1994); C. Tong and Z. Warhaft, Phys. Fluids 6, 2165 (1994). 
[5] R. H. Kraichnan, Phys. Fluids 11, 945 (1968); J. Fluid. Mech. 64, 737 (1974); Phys. 

Rev. Lett. 72, 1016 (1994); ibid. 78, 4922 (1997). 
[6] M. Chertkov, G. Falkovich, I. Kolokolov and V. Lebedev, Phys. Rev. E 52, 4924 (1995); 

M. Chertkov and G. Falkovich, Phys. Rev. Lett. 76, 2706 (1996). 
[7] K. Gawedzki and A. Kupiainen, Phys. Rev. Lett. 75, 3834 (1995); D. Bernard, K. 

Gawedzki and A. Kupiainen, Phys. Rev. E 54, 2564 (1996). 
[8] B. I. Shraiman and E. D. Siggia, Phys. Rev. Lett. 77, 2463 (1996); A. Pumir, B. I. 

Shraiman and E. D. Siggia, Phys. Rev. E 55, R1263 (1997); A. Pumir, Europhys. Lett. 

37, 529 (1997); Phys. Rev. E 57, 2914 (1998). 
[9] L. Ts. Adzhemyan, N. V. Antonov and A. N. Vasil'ev, Phys. Rev. E 58, 1823 (1998). 
[10] M. Chertkov, G. Falkovich, and V. Lebedev, Phys. Rev. Lett. 76, 3707 (1996). 
[11] G. Eyink, Phys. Rev. E 54, 1497 (1996). 

[12] N. V. Antonov, e-Print archive LANL, \zhao-dyn/9808011{ , submitted to Phys. Rev. E. 
[13] T. Elperin, N. Kleeorin, and I. Rogachevskii, Phys. Rev E 52, 2617 (1995); ibid. 55, 

2713 (1997); Phys. Rev. Lett. 76, 224 (1996). 
[14] M. Vergassola and M. Avellaneda, Physica D 106, 148 (1997). 
[15] K. Gaw§dzki and M. Vergassola, e-Print archive LANL, \cond-mat/9811 3% . 
[16] A. Celani, A. Lanotte and A. Mazzino, e-Print archive LANL, \chao-dyn/99U20Tl . 



[17] L. Ts. Adzhemyan and N. V. Antonov, Phys. Rev. E 58, 7381 (1998) 

[18] M. Vergassola and A. Mazzino, Phys. Rev. Lett. 79, 1849 (1997). 

[19] A. Lanotte and A.Mazzino, e-Print archive LANL, \chao-dyn/9903026 . 

[20] I. Arad, B. Dhruva, S. Kurien, V. S. L'vov, I. Procaccia and K. R. Sreenivasan, Phys. 
Rev. Lett, 81, 5330 (1998); I. Arad, L. Biferale, I. Mazzitelli and I. Procaccia, e-Print 
archive LANL, \chao-dyn/9903013i , S. Kurien, V. S. L 'vov, I. Procaccia, and K.R. Sreeni- 
vasan, e-Print archive LANL, \;hao-dyn/990603~l . 

[21] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasil'ev, Usp. Fiz. Nauk, 166, 1257 (1996) 
[Phys. Usp. 39, 1193 (1996)]. 

[22] L. Ts. Adzhemyan, N. V. Antonov and A. N. Vasiliev, The Field Theoretic Renormal- 
ization Group in Fully Developed Turbulence (Gordon & Breach, London, 1999). 

[23] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Clarendon, Oxford, 
1989). 

[24] J. P. Bouchaud, A. Comtet, A. Georges, and P. Le Doussal, J. Phys. (Paris) 48, 1445 

(1987); 49, 369 (1988). 
[25] J. P. Bouchaud and A. Georges, Phys. Rep. 195, 127 (1990). 

[26] J. Honkonen and E. Karjalainen, J. Phys. A: Math. Gen. 21, 4217 (1988); J. Honko- 



23 



nen, Yu. M. Pis'mak, and A. N. Vasil'ev, J. Phys. A: Math. Gen. 21, L835 (1989); J. 
Honkonen and Yu. M. Pis'mak, J. Phys. A: Math. Gen. 22, L899 (1989). 
[27] L. Ts. Adzhemyan, S. V. Borisenok, and V. I. Girina, Theor. Math. Phys. 105, 1556 
(1995). 

[28] L. Ts. Adzhemyan, N. V. Antonov, M. Hnatich and S. V. Novikov, 'Anomalous scaling 
of a passive scalar in the presence of strong anisotropy," in preparation. 



24 



TABLES 



TABLE I. Canonical dimensions of the fields and parameters in the model (2.7). 
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